Large-scale binding affinity calculations on commodity compute clouds

In recent years, it has become possible to calculate binding affinities of compounds bound to proteins via rapid, accurate, precise and reproducible free energy calculations. This is imperative in drug discovery as well as personalized medicine. This approach is based on molecular dynamics (MD) simulations and draws on sequence and structural information of the protein and compound concerned. Free energies are determined by ensemble averages of many MD replicas, each of which requires hundreds of cores and/or GPU accelerators, which are now available on commodity cloud computing platforms; there are also requirements for initial model building and subsequent data analysis stages. To automate the process, we have developed a workflow known as the binding affinity calculator. In this paper, we focus on the software infrastructure and interfaces that we have developed to automate the overall workflow and execute it on commodity cloud platforms, in order to reliably predict their binding affinities on time scales relevant to the domains of application, and illustrate its application to two free energy methods

Combining mental imagery and implementation intentions to increase brisk walking in people with the symptoms of depression

Volume one consists of three parts. Part one is a meta-analytic review investigating the effect of various doses of physical exercise on the symptoms of depression. The quality of twelve Randomised Controlled Trials (RCTs) was assessed. Outcomes from four RCTs were included in a meta-analysis. Results indicated that higher doses of exercise were more effective than no exercise in reducing the symptoms of depression, whereas there was no difference between lower doses of exercise and no exercise. The findings are interpreted with caution and more robust research is required to investigate the relationship between exercise and depression, addressing the methodological limitations discussed in this review. Part two is an empirical study investigating the extent to which people with the symptoms of depression could generate mental imagery to increase goal attainment. Sixty five students were given the goal of increasing brisk walking over two weeks. According to randomisation participants were either asked to rehearse the goal, to form implementation intentions about the goal or to imagine themselves undertaking the goal. Brisk walking increased with no between-group differences. It was tentatively concluded that being given a goal intention was sufficient to improve attainment. Further research, controlling for research participation effects and natural changes in behaviour, is required to confirm the finding. Part three is a critical appraisal of the meta-analytic review and empirical study. It explores the personal implications the empirical study has had upon clinical work. It also evaluates the empirical study’s web-based design, highlighting the methodological and ethical issues that were raised. Finally, the process of undertaking a meta-analysis is discussed, with specific focus on the challenges posed by this methodology

Atomistic Modeling of Scattering Curves for Human IgG1/4 Reveals New Structure-Function Insights

Small angle x-ray and neutron scattering are techniques that give solution structures for large macromolecules. The creation of physically realistic atomistic models from known high-resolution structures to determine joint x-ray and neutron scattering best-fit structures offers a, to our knowledge, new method that significantly enhances the utility of scattering. To validate this approach, we determined scattering curves for two human antibody subclasses, immunoglobulin G (IgG) 1 and IgG4, on five different x-ray and neutron instruments to show that these were reproducible, then we modeled these by Monte Carlo simulations. The two antibodies have different hinge lengths that connect their antigen-binding Fab and effector-binding Fc regions. Starting from 231,492 and 190,437 acceptable conformations for IgG1 and IgG4, respectively, joint x-ray and neutron scattering curve fits gave low goodness-of-fit R factors for 28 IgG1 and 2748 IgG4 structures that satisfied the disulphide connectivity in their hinges. These joint best-fit structures showed that the best-fit IgG1 models had a greater separation between the centers of their Fab regions than those for IgG4, in agreement with their hinge lengths of 15 and 12 residues, respectively. The resulting asymmetric IgG1 solution structures resembled its crystal structure. Both symmetric and asymmetric solution structures were determined for IgG4. Docking simulations with our best-fit IgG4 structures showed greater steric clashes with its receptor to explain its weaker FcγRI receptor binding compared to our best-fit IgG1 structures with fewer clashes and stronger receptor binding. Compared to earlier approaches for fitting molecular antibody structures by solution scattering, we conclude that this joint fit approach based on x-ray and neutron scattering data, combined with Monte Carlo simulations, significantly improved our understanding of antibody solution structures. The atomistic nature of the output extended our understanding of known functional differences in Fc receptor binding between IgG1 and IgG4

Domain structure of human complement C4b extends with increasing NaCl concentration: implications for its regulatory mechanism

During the activation of complement C4 to C4b, the exposure of its thioester domain (TED) is crucial for the attachment of C4b to activator surfaces. In the C4b crystal structure, TED forms an Arg(104)-Glu(1032) salt bridge to tether its neighbouring macroglobulin (MG1) domain. Here, we examined the C4b domain structure to test whether this salt bridge affects its conformation. Dual polarisation interferometry of C4b immobilised at a sensor surface showed that the maximum thickness of C4b increased by 0.46 nm with increase in NaCl concentration from 50 mM to 175 mM NaCl. Analytical ultracentrifugation showed that the sedimentation coefficient s20, w of monomeric C4b of 8.41 S in 50 mM NaCl buffer decreased to 7.98 S in 137 mM NaCl buffer, indicating that C4b became more extended. Small angle X-ray scattering reported similar RG values of 4.89-4.90 nm for C4b in 137-250 mM NaCl. Atomistic scattering modelling of the C4b conformation showed that TED and the MG1 domain were separated by 4.7 nm in 137-250 mM NaCl, this being greater than that of 4.0 nm in the C4b crystal structure. Our data reveal that in low NaCl concentrations, both at surfaces and in solution, C4b forms compact TED-MG1 structures. In solution, physiologically-relevant NaCl concentrations lead to the separation of the TED and MG1 domain, making C4b less able to bind to its complement regulators. These conformational changes are similar to those seen previously for complement C3b, confirming the importance of this salt bridge for regulating both C4b and C3b

An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics

The accurate identification of the specific points of interaction between G protein-coupled receptor (GPCR) oligomers is essential for the design of receptor ligands targeting oligomeric receptor targets. A coarse-grained molecular dynamics computer simulation approach would provide a compelling means of identifying these specific protein–protein interactions and could be applied both for known oligomers of interest and as a high-throughput screen to identify novel oligomeric targets. However, to be effective, this in silico modeling must provide accurate, precise, and reproducible information. This has been achieved recently in numerous biological systems using an ensemble-based all-atom molecular dynamics approach. In this study, we describe an equivalent methodology for ensemble-based coarse-grained simulations. We report the performance of this method when applied to four different GPCRs known to oligomerize using error analysis to determine the ensemble size and individual replica simulation time required. Our measurements of distance between residues shown to be involved in oligomerization of the fifth transmembrane domain from the adenosine A2A receptor are in very good agreement with the existing biophysical data and provide information about the nature of the contact interface that cannot be determined experimentally. Calculations of distance between rhodopsin, CXCR4, and β1AR transmembrane domains reported to form contact points in homodimers correlate well with the corresponding measurements obtained from experimental structural data, providing an ability to predict contact interfaces computationally. Interestingly, error analysis enables identification of noninteracting regions. Our results confirm that GPCR interactions can be reliably predicted using this novel methodology

Laser Generation and Detection of Surface Acoustic Waves Using Gas-Coupled Laser Acoustic Detection

Laser generation and detection of ultrasound has the advantage of requiring no mechanical contact with the materials under investigation. We previously reported [1] laser-based measurements on Lamb waves in graphite/polymer composite laminates using a confocal Fabry-Perot interferometer for detection. Related work by other groups includes air-coupled detection of Lamb waves in similar composites using capacitive transducers [2,3] and interferometric detection of Lamb waves in paper [4]. Our earlier work has been extended using Gas-Coupled Laser Acoustic Detection (GCLAD), an economical alternative laser-based method which has the additional advantage that the detection laser beam is not reflected from the sample surface. GCLAD is thus particularly useful for materials with surfaces of poor optical quality. We demonstrate below that the combination of laser generation and GCLAD can be used to obtain well-resolved surface-acoustic waves (SAWs) in a variety of materials, including metals, paper, thin films, and composite pre-preg tape. We also show some preliminary SAW scans obtained with laser generation and GCLAD using metallic samples. Each pixel in the scans represents the strength of a SAW passing through a portion of the sample with an area of about 1 cm2. Scans of this type offer the possibility of economical testing of large sample areas, potentially on-line in a manufacturing environment

Application of the ESMACS Binding Free Energy Protocol to a Multi‐Binding Site Lactate Dehydogenase A Ligand Dataset

Over the past two decades, the use of fragment‐based lead generation has become a common, mature approach to identify tractable starting points in chemical space for the drug discovery process. This approach naturally involves the study of the binding properties of highly heterogeneous ligands. Such datasets challenge computational techniques to provide comparable binding free energy estimates from different binding modes. The performance of a range of statistically robust ensemble‐based binding free energy calculation protocols, called ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent), is evaluated. Ligands designed to target two binding pockets in the lactate dehydogenase, a target protein, which vary in size, charge, and binding mode, are studied. When compared to experimental results, excellent statistical rankings are obtained across this highly diverse set of ligands. In addition, three approaches to account for entropic contributions are investigated: 1) normal mode analysis, 2) weighted solvent accessible surface area (WSAS), and 3) variational entropy. Normal mode analysis and WSAS correlate strongly with each other—although the latter is computationally far cheaper—but do not improve rankings. Variational entropy corrects exaggerated discrimination of ligands bound in different pockets but creates three outliers which reduce the quality of the overall ranking

Non-linearity of the collagen triple helix in solution and implications for collagen function

Collagen adopts a characteristic supercoiled triple helical conformation which requires a repeating (Xaa-Yaa-Gly)n sequence. Despite the abundance of collagen, a combined experimental and atomistic modelling approach has not so far quantitated the degree of flexibility seen experimentally in the solution structures of collagen triple helices. To address this question, we report an experimental study on the flexibility of varying lengths of collagen triple helical peptides, composed of six, eight, ten and twelve repeats of the most stable Pro-Hyp-Gly (POG) units. In addition, one unblocked peptide, (POG)10unblocked, was compared with the blocked (POG)10 as a control for the significance of end effects. Complementary analytical ultracentrifugation and synchrotron small angle X-ray scattering data showed that the conformations of the longer triple helical peptides were not well explained by a linear structure derived from crystallography. To interpret these data, molecular dynamics simulations were used to generate 50 000 physically realistic collagen structures for each of the helices. These structures were fitted against their respective scattering data to reveal the best fitting structures from this large ensemble of possible helix structures. This curve fitting confirmed a small degree of non-linearity to exist in these best fit triple helices, with the degree of bending approximated as 4–17° from linearity. Our results open the way for further studies of other collagen triple helices with different sequences and stabilities in order to clarify the role of molecular rigidity and flexibility in collagen extracellular and immune function and disease

Radial Velocities as an Exoplanet Discovery Method

The precise radial velocity technique is a cornerstone of exoplanetary astronomy. Astronomers measure Doppler shifts in the star's spectral features, which track the line-of/sight gravitational accelerations of a star caused by the planets orbiting it. The method has its roots in binary star astronomy, and exoplanet detection represents the low-companion-mass limit of that application. This limit requires control of several effects of much greater magnitude than the signal sought: the motion of the telescope must be subtracted, the instrument must be calibrated, and spurious Doppler shifts "jitter" must be mitigated or corrected. Two primary forms of instrumental calibration are the stable spectrograph and absorption cell methods, the former being the path taken for the next generation of spectrographs. Spurious, apparent Doppler shifts due to non-center-of-mass motion (jitter) can be the result of stellar magnetic activity or photospheric motions and granulation. Several avoidance, mitigation, and correction strategies exist, including careful analysis of line shapes and radial velocity wavelength dependence.Comment: Invited review chapter. 13pp. v2 includes corrections to Eqs 3-6, updated references, and minor edit